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SHOCK TUBE MEASUREMENTS OF SPECIFIC REACTION 
RATES IN THE BRANCHED CHAIN CH 4 -C0-0 2 SYSTEM 

T. A. Brabbs and R. S. Brokaw 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

Rate constants of two elementary bimolecular reactions involved 
in the oxidation of methane have been determined by monitoring the 
exponential growth of CO flame band emission behind incident shocks 
in three suitably chosen gas mixtures. Results are as follows: 

O + CH 4 -CH 3 +OH, k 7 = 1.9 x 10 14 exp (-5900/T) cm 3 mole" 1 sec" 1 
(1300 - 2000K) andCHg + 0 2 ~CHgO +0, k 1Q = 2. 4 x 10 13 exp(-14500/T) 
(1200 - 1800K) assuming CHgO + M-CH 2 0 + H + M. (If it is assumed 
that CHgO + CO-C0 2 + CH 3 , then k 1Q = 5.9 x 10 13 exp(-14500/T).) 

The data do not support a mechanism which invokes the four center 
process CHg + Og-CHgO + OH for the reaction of methyl with oxygen. 


INTRODUCTION 


When a mixture containing carbon monoxide, oxygen, and a small 
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amount of methane is subject to a temperature and pressure pulse 
in a shock tube, small concentrations of atoms and free radicals 
are first formed by processes such as 

CO + 0 2 —C0 2 + 0« 

These atom and radical concentrations may then grow exponentially 
via the branched chain scheme 

k 0 

h + o 2 -OH + O (II) 

ke 

OH + CO - C0 2 + H (V) 

k„ 

o+ch 4 -^ch 3 + oh 0 (vn) 

Chains may be terminated by the reactions 

k 

H + 0 2 + M 4. H0 2 + M (IV) 

k fi 

O + CO + M 2. C0 2 + M 0 (VI) 

The reactions 

k 8 

H + ch 4 a 

k 9 

OH+CH 4 a 


h 2 + ch 3 (vm) 

h 2 o + ch 3 


(IX) 
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may also serve as chain terminations under circumstances where 
methyl radicals are unreactive . (There reactions have been 
numbered to be compatible with notation used in a previous study 
of the H 2 - CO - Og system. *) 

Theory shows 2, 3 that in such chain-branched systems the 
atom and radical concentrations increase as expXt (except very 
early in the reaction) . The growth constant X depends on the rate 
constants of the elementary chemical reactions and the concentrations 
of stable reactants- -in this case CO, 0 2 , and CH^. 

This paper reports experimental growth constants measured behind 
incident shocks for a range of pressures, temperatures and gas 
compositions. It is found that reactions (n) - (IX) do not suffice to 
explain the experimental results; rather, it is necessary to consider 
reaction of the methyl radical with oxygen. The data are consistent 
with reactions (II) - (IX) and with methyl reacting according to 

ch 3 + o 2 kl 9- ch 3 o + O (X) 


followed by either 


CH 3 0 + M X i- CH 2 0 + H + M 


or 

^1 o 

CHgO + CO CH 3 + C0 2 . 

The data are not consistent with the four center reaction 

If* 

“10 xt r\ . 


(XI) 

(xn) 

t 

on 


ch 3+ ° 2 


CH 2 0 + OH 
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which has often been proposed 4 for the reaction of methyl radicals 
in combustion processes, 

- EXPERIMENTAL ASPECTS 

Growth constants were obtained by measuring the blue CO flame 
band emission behind incident shocks. The intensity of this radiation 
is proportional to the product of carbon monoxide and oxygen atom 
concentrations, 5 and since very little CO is consumed, the light 
measures the increase of oxygen atom concentration with time. 

Details of the shock tube and associated optical and electronic equip- 

g 

ment have been described elsewhere. 

Gas mixtures contained small amounts of CH^ with varying amounts 
of C0,0 2> and in one mixture COg, diluted with argon. (Carbon dioxide 
was added to ensure vibrational relaxation of carbon monoxide in the 
mixture with a large CO concentration.) Methane, oxygen, and argon 
were high-purity tank gases and were used without further purification. 
Carbon monoxide and carbon dioxide were purified as previously 
described. * 

fi 7 

It has been shown ’ that boundary layer effects must be considered 
in analyzing data obtained behind incident shocks. Conditions behind 
the shocks, in the region of the experimental measurements, were 
obtained from a computer program which integrates the equations of 
chemical change for a shocked gas accounting for the effects of boundary 
layer buildup; the procedure has been described in reference 7 for 
laminar boundary layers. In this work the boundary layer was turbulent; 
the maximum separation distance between shock and contact surface 
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was calculated according to reference 9. 

Exponential growth constants were obtained from plots of the 
logarithm of observed light intensity versus gas time. The 
relation between gas and laboratory times was obtained from the 
computer calculations. 

THEORETICAL CONSIDERATIONS 

The analytic solution of the differential equations describing 

the ignition kinetics in chain-branched systems involving hydrogen, 

oxygen and carbon monoxide has been discussed in detail in previous 
2 3 11 

publications. 9 * The solution for the system of equations 
discussed in the introduction is quite analogous and will merely be 
outlined here. 

To obtain the analytic solution the atom and radical concentrations 
are assured to be small in comparison to the concentrations of 
reactants so that the reactant concentrations may be considered 
constant and reactions between chain carriers are unimportant. In 
addition the temperature and pressure must be nearly constant over 
the time range of the experimental observations. That these assumptions 

O 

were in fact satisfied was confirmed by numerical computations. 

Subject to these assumptions the chemical kinetics are described 
by a system of simultaneous first order linear differential equations, 
one for each chain carrier. The solution of these equations shows that 
the chain carrier concentrations grow as expAt (except very early in 
the reaction). The growth constant A is the positive root of a polynomial 

of a degree dependent on the number of chain carriers. Thus with 
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H, OH, O, CHg, and CH^O as chain carriers a quintic equation 
is obtained. Actually there is little information in the literature 
as to rates of reaction of the methoxy radical. We carried out 
numerical calculations assuming reasonable but rapid rates for 
reactions (XT) and (XU). (The rate constant for reaction (XI) was 
assumed to be kjj = 5 x 10*^ exp (-10, 500/T), based on a reasonable 
preexponential factor with an activation energy set equal to the 
endothermicity of the reaction. The rate constant for reaction (XU) 
was assumed equal to the known rate for the analogous reaction 
(V)) The numerical calculations showed that the relation between 
CHgO and CHg concentrations is in reasonable agreement with the 
steady state approximation, 

k 10[°2H CH 3] = (knM +k 12 [C0])[CH 3 O] (1) 

Equation (1) was substituted into the differential equations describing 
the chemical kinetic to eliminate the methoxy radical concentration. 
This reduced the polynomial in X to the quartic equation 

-(v 2 + v^+v q +X) u 5 0 

v 2 -(^ + j, 9+ X) 

0 -(i/g + i^+X) 

"8 "9 “7 

Where 1^5 1^[0 2 ], u^S k 4 [0g] [M] , i^®k 5 [CO], 
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^5 kg [CO] [M], ^k 7 [CH 4 ], ^g^kgtCH^, 

^ sk g[CH 4 ], ^0 = k io[°23> ^10 " k lo[°2^ and 
aSkn[M] / (k^[M] + kj 2 [CO]) 0 If the methoxy radical reacts 

only via reaction (XI) ar = lj if it reacts only via reaction (XII), 

a = 0 . 

Equation (2) was used in the data analysis to obtain rate constants 
for the individual reactions k 7 , k jq, and k 1Q . For example, it can be 
solved for Uj and hence k 7 to obtain an expression involving the rates 
of the other reactions in the scheme and the experimental growth 
constant. 


RESULTS AND DISCUSSION 

The compositions of the three gas mixtures which were investigated 
are given in Table I . (Table I also has calculated sensitivities which 
show the importance of each reaction in determining the growth constant.) 
Initially two mixtures, A and B, were studied. In mixture A the 
concentrations of methane and oxygen are small, so that the growth 
constant should depend principally on the rates of reactions (II) and 
(VII). On the other hand in mixture B the methane and carbon monoxide 
concentrations are small so that the growth constant should depend 
principally on the rates of reactions (V) and (VH). This is confirmed 
by the growth constant sensitivities shown in Table I. 

Initially it was assumed that the methyl radical did not react with 
either carbon monoxide or oxygen. In this event the expression for 
the growth constant reduces to a cubic, the determinant of the first 
three rows and columns of equation (2) equated to zero. The cubic was 
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solved for k^ and values of k ? were calculated from experimental 
growth constants for mixtures A and B, with rate constants for 
the other reactions shown in Table II. The rate constants for 
reactions (II) and (V), which are important, were determined in 
our apparatus. Rate constants for reactions (IV), (VI), (VHI) and 
(IX) have less effect on growth constants; they are in the nature of 
corrections. Reactions (IV) and (VI) are completely negligible 
except at the lower temperatures. The rate for reaction (IV) was 
deduced from unpublished shock tube experiments at Lewis Research 
Center by M. R. Lauver. This value, obtained near HOOK, is 
consistent with the data obtained in this paper inasmuch as the rate 
constants of reference 1 were used in obtaining k^. (Third body 
efficiencies were taken from reference 14.) 

Rate constants computed for reaction (VII) are shown in figure 1. 

Values computed from mixture A are in rough agreement with an 

15 

extrapolation of the low temperature data of Wong and Potter. 
However, rate constants derived from mixture B are substantially 
larger, particularly at the higher temperatures. Mixture B contained 
a large amount of oxygen; thus it appears that the methyl radical is 
in fact reacting with oxygen and continuing the chain. Consequently 
we expanded our analysis to account for reaction of the methyl radical 
according to reactions (X*) and (X); the methoxy radical produced in 
reaction (X) may either be decomposed in reaction (XI) ( a = 1) or 
or react with carbon monoxide in reaction (XII) (a = 0). 

In order to investigate these possible reaction paths mixture C was 
prepared, similar to mixture B except for a tenfold increase in 



9 


methane concentration. Equation (2) was solved for k^Q and kjQ 
and the resulting expressions were used to compute these rate 
constants from the experimental growth constants using other rate 
constants cited in table II along with a rate for reaction (VII) from 
reference 15. 

Rate constants derived from mixture B and C assuming methyl 
and oxygen react via the four center reaction (X") are shown in 
figure 2(a) . Rate constants from the two mixtures differ by almost 
an order of magnitude. Further if the Arrhenius equation is 
applied to the less scattered data from mixture C, the preexponential 
factor is about 2 X 10 17 — about three orders of magnitude in excess 
of collision frequency! Thus we conclude the reaction of methyl 
with oxygen via reaction (X"*) cannot describe our results. 

Rate constants derived assuming methyl and oxygen react via 
reaction (X) to form methoxy and atomic oxygen are shown in 
figure 2(b) and (c) . The data scatter much less than in figure 2(a) . 
Furthermore the preexponential Arrhenius factor is reasonable; 

2.4X 10* 3 if CHgO decomposes via reaction (XI) (a = 1), or 
5. 9 x 10 13 if it reacts with carbon monoxide via reaction (XII) 

(a = 0). Finally, the activation temperature E/R, about 14500K, 
is very close to a minimum value of 14600K based on the enjothermicity 
of reaction (X) 16 . Thus it seems probable that methyl and oxygen 
react via the simple abstraction process (X), although it is not 
possible to reach a conclusion from figure 2(b) and (c) as to the fate 
of methoxy. 
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In an attempt to resolve this point, rate constants for reaction (X) 

(the lines of figure 2(b) and (c)) derived from mixture C were used 

to calculate rates for reaction (V) from mixture B which were then 

compared with values previously obtained from a study of the 

H 2 -C0-0 2 system. 1 Agreement was slightly better for values 

calculated assuming decomposition of methoxy, reaction (XI), 

especially as to the trend with temperature. Furthermore, there 

is evidence that at room temperature the reaction between carbon 

17 -3 “4 

monoxide and methoxy is slow, with a rate 10 to 10 of 

18 

reaction (V) between carbon monoxide and hydroxyl. Thus we 
favor reaction (XI) over reaction (XII) and have assumed this path 
in our final data analysis, although reaction (XII) can also adequately 
account for our data. 

Table m presents experimental temperatures and pressures 
(corrected for boundary layer effects) with growth constants for 
mixtures A, B, and C. Also shown are rate constants for reactions 
(VII) and (X) based on the final data analysis. These values were 
obtained as follows: Equation (2) was solved to obtain expressions 
for k ? and k 1Q , with k' Q = 0 and a - 1. A trial value of k 1Q from 
figure 2(b) was used with other rate constants from Table II to compute 
values of k^ from growth constants from mixtures A and B. The 
resulting rate constants were fitted to an Arrhenius equation and used 
with growth constants from mixture C to compute values of k^g. 

These values of k^g were in turn fitted to an Arrhenius equation and 
used to recompute values of k ? . The process was repeated until 
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values of and k 1Q no longer changed. (One recomputation of 
k^ andkjQ sufficed.) 

In contrast to the results shown in figure 1, values of k^ from 
mixtures A and B were similar, and are shown in figure 3. A least 
squares fit to the data for both mixtures gives 

k ? = 1,9 X 10 14 exp(-5900/T) cm 3 mole -1 sec' 1 (3) 

in the temperature range 1300 - 2000K, in reasonable agreement with 
the low temperature data of reference 15. The standard deviation in 
the activation temperature is 570K and the standard deviation in lnk^ 
is 0.24, or about 28%. 

A least squares fit to the rate constants for reaction (X) from 
mixture C is also shown in figure 3. The rate constant is 

kjQ = 2. 4 x 10 13 exp (-14500/T) cm 3 mole -1 sec 1 (4) 

in the temperature range 1200 - 1800K. The standard deviation in 
the activation temperature is 870K and the standard deviation in lnkj Q 

is 0.15, or about 16%. Equation (4) was derived assuming destruction 
of methoxy via reaction (XI) ; if methoxy is consumed by reaction (XU) 

13 

the preexponential factor in equation (4) should be increased to 5.9 x 10 . 

Finally the sensitivities of the growth constants to the various rate 
constants were computed numerically from equation (2) and tabulated 
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in Table I. In general, growth constants depend most strongly on 
the rates of reactions (II), (V), (VII) and (X) which have been 
determined in the shock tube used here. However, there is 
appreciable sensitivity to the rates of reactions (Vm) and (IX), 
taken from reference 13. If in the future high temperature studies 
indicate appreciably different rates for reactions (VTII) and (IX) 
a reanalysis of the data obtained in this paper would be in order 
to obtain improved values of k ? and k 1Q . We very much doubt 
that any changes in kg and kg could alter our conclusion as to the 
course of the reaction between the methyl radical and oxygen. 

CONCLUDING REMARKS 

An important result of this paper is the conclusion that in the 
high temperature oxidation of methane the methyl radical reacts with 
oxygen to form methoxy and atomic oxygen, rather than formaldehyde 
and hydroxyl as has frequently been assumed heretofore. 

In addition, high temperature reaction rates have been obtained 
for the reactions O + CH 4 ~ CHg + OH and CHg + 0 2 - CHgO + O, in 
the range 1200 - 2000K, 
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TABLE I. 

MIXTURE COMPOSITIONS AND GROWTH CONSTANT SENSITIVITIES 

Mixture 

ABC 


Composition, % 

ch 4 

0.0194 

0.0193 

0.205 

CO 

18.72 

1.032 

1.056 

°2 

1.033 

19.41 

19.54 

co 2 

5.68 

- 

- 

Ar 

74.5 

79.5 

79.2 

Sensitivities (150 OK) 

BlnX/Slni^ 

0.74 

0.06 

0.07 

01nX/31m^ 

-.07 

-.03 

-.02 

01nX/01ni^ 

.12 

.47 

. 30 

01nX/01n^g 

-.12 

.00 

.00 

01nX/01ni^ 

.56 

.41 

.15 

81nX/01m^ 

-.19 

o 

o 

• 

-.03 

01n X/01mjj 

-.05 

-.15 

-. 20 

01nX/31ni^Q 

.02 

.24 

.74 
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TABLE H. 

RATE CONSTANTS USED IN DATA ANALYSIS 



Reaction 

Rate Constant 

(mole, cm, sec, K units) 

Reference 

(II) 

h+o 2 -oh+o 

1.38xiO l4 exp(-8250/T) 

7 

(IV) 

H +0 2 + M — HOg + M 

10 15 (M=Ar) 

See Text 

(V) 

oh + co-co 2 +h 

4. 2X10 11 exp(-500/T) 

1 

(VI) 

o+co + m-co 2 + m 

3x 10 14 exp (-1510/T) 

12 

(vn) 

H + CH 4 -H 2 + CH 3 

1. 24x 10 14 exp(-6000/T) 

13 

(IX) 

oh+ch 4 -h 2 o+ch 3 

2„ 5x 10 13 exp(-2520/T) 

13 
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TABLE m. - EXPERIMENTAL RESULTS 


T(K) 

P(atm) 

X (sec “ 3 ) 

3 - i 

k (qm mole' 

Mixture A 
1200 

1.26 

3.01 x 10 3 

k 7 

2. 14 x 10 3 

1302 

1.39 

4.92 

2.08 

1372 

.99 

6. 89 

3.41 

1498 

1.11 

9.84 

3.54 

1519 

1.12 

10. 56 

3.69 

1546 

1.15 

11.46 

3.84 

1672 

.94 

14. 08 

. 5.93 

1719 

.97 

15.45 

6.17 

1817 

.87 

17.99 

8.91 

1990 

.78 

23.03 

15. 63 

Mixture B 
1221 

1.84 

5.71 x 10 3 

k ? 

1.15 x 10 1 

1255 

1,69 

7.09 

1.63 

1281 

1.62 

8.25 

2.13 

1347 

1.49 

8,93 

2.36 

1353 

1.62 

8.59 

1.85 

1404 

1.45 

9.80 

2.60 

1434 

1.37 

9.29 

2.41 

1444 

1.00 

7.57 

2.91 

1460 

1.26 

12.51 

5.49 

1515 

1.06 

10. 42 

4,57 

1515 

1.06 

9.97 

4. 09 

1536 

1.21 

11.02 

3.59 

1578 

1.11 

10.01 

3. 14 

1624 

1.15 

11.40 

3.53 

1670 

1.04 

12.18 

4.76 

1762 

.95 

12.51 

5.19 

Mixture C 
1278 

2.00 

8.11 X 10 3 

k 10 

. 37 X 10 9 

1407 

1.48 

9.80 

.68 

1431 

1.75 

12.58 

,77 

1464 

1.42 

12.31 

.99 

1499 

1.32 

15.25 

1.46 

1499 

1.59 

17.92 

1.43 

1546 

1.41 

24. 24 

2.11 

1587 

1.50 

22.36 

2.61 

1601 

1.14 

20.20 

2.73 

1654 

1.33 

31.12 

4.16 

1677 

1.20 

29.90 

4. 58 
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Figure 1. - Rate constants for reaction (VII) assuming methyl unreactive 
with CO and 0 2 . 



<a|CH 3 + 0 2 - CH^ +0K 


TEMPERATURE, K 

(b) CH 3 + 0 2 - CH 3 O + 0; CH 3 O + M - CH 2 0 + H + M. (c) CH 3 + 0 2 - CH 3 0 + 0; CH 3 O + CO - C0 2 + ch 3 

Figure 2 . - Rate constants for reaction of CH 3 with 0 2 . 
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